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1. Introduction and summary 

A black hole in an isolated cavity which is bigger than its Schwarzschild radius can be in equilibrium 
with surrounding radiation only if the total energy E of the system is greater than a critical value Eq, 
which depends on the volume of the cavity and the number of fields in the radiation (Hawking 1976). 
While local equilibrium configurations exist for energies greater than Ec, it is generally assumed that 
black holes must evaporate if Eq < E < Eb — l-3£'c since for E < Eb pure radiation has more 
entropy than the composite system (Gibbons and Perry 1978). 

This holds only if the system is given unlimited time to relax. If one considers instead a sequence of 
configurations starting, say, from an equilibrium at E = Eb and in which the total energy is decreased 
quasi-statically, the system reaches new local equilibrium configurations in a finite time (Zurek 1980). 
During that time, the probability of finding a fluctuation high enough that the black hole completely 
evaporates is exponentially small. So small that the black hole is most likely to survive in a metastable 
superheated state up to energies very close indeed to Ec (Okamoto, Katz and Parentani 1994). This 
is true for cavities whose radius L* is greater than 10^ Planck lengths Ip = {fiG / c^Y^'^ ~ 10^^'^cm 
(L = L* /Ip > 10^). For those cavities, back-reaction due to quantum matter effects (York 1985), 
quantum gravity and spontaneous nucleation due to thermal fluctuations (Piran and Wald 1982) 
may be neglected since the maximum temperature, reached by equilibrium conflgurations, is always 
smaller than To ~ 0.37-L~^/^Q (in Planck units). Thus more than three orders of magnitude separate 
the equilibrium temperature from the Planck one. This indicates that a thermodynamic analysis may 
be safely performed. On the contrary, for smaller boxes the validity of this analysis becomes dubious. 

Self-gravitating thermal radiation exhibits a similar behavior in an evolutionary sequence in which 

the total energy increases quasi-statically. Indeed, there is now an upper limit of energy Ea — 0.25L 

(here E is measured in Planck units) above which no equilibrium conflguration exists (Klein 1947). 

Thus for energies Eb < E < Ea {Eb ~ 0.20^3/5 <^ Ea ^ 0.25L for L > 10^) pure radiation is in 

metastable superheated states. We shall prove that, as in the black hole situation, pure radiation is 

most likely to remain in those metastable states almost up to Ea- At that point, a thermodynamic 

as well as dynamical instability (Sorkin, Wald and Zhang 1980) develops and most of the radiation 
* Throughout the text numerical values are given with two significant digits 



collapses into a black hole leading to a composite state of a black hole in equilibrium with surrounding 
left-over radiation. The spontaneous evaporation of the black hole into pure radiation near Ec and 
the collapse of the radiation near Ea provides thus a unified picture relating the different phases of a 
black hole and self-gravitating radiation in a cavity. 

Uniform rotation added to the cavity gives a total angular momentum J shared between the 
radiation and the black hole. The equilibrium configurations change in the following way. The critical 
energy Ec{J) at which the black hole evaporates is now a function of J. It increases with increasing 
angular momentum while the fraction of black hole energy at that point E},hc/Ec decreases from 4/5 
at J = to almost zero in the limiting rotating case: J = Hence for fast rotating cavities, 

at fixed J and with slowly decreasing total energy, superheated Kerr black holes will evaporate with 
little latent heat, since black holes survive almost down to the critical point Ec{J)- 

All of this concerns microcanonical situations in which the total energy is the control parameter. 
In canonical situations, with the temperature kept fixed, almost all the configurations wherein a black 
hole coexists with radiation are unstable. There exits nevertheless a narrow range of energies for very 
fast rotating cavities in which the canonical ensemble is stable. One thus recovers in this narrow range 
the already noticed flip of the heat capacity for Kerr black holes (Davies 1981, Kaburaki et al. 1993; 
cf. Katz et al. 1993 for Kerr- Newman black holes). When the radiation is included, this happens, 
however, when gravitational effects (neglected in the present analysis) are important. 

In this paper, we analyze the stability limits, the different degrees of instability and the fluctuations 
of each of the two states of the system (with or without black hole) in evolutionary sequences of equi- 
librium configurations, with and without angular momentum. The analysis is purely thermodynamic. 
We assume there exists a state function whose stationary value is a function of the total mass-energy 
E, total angular momentum J and the volume of the cavity V = XL^. We consider separately the mi- 
crocanonical and the canonical ensembles since they are not equivalent. In a microcanonical ensemble, 
E is the main control parameter and the inverse temperature 

kT dE ' ' ' ^ ' 

is the conjugate parameter of E with respect to the total entropy S. In a canonical ensemble, the 
Massieu function is not S but S/k — (3E = —fiF (F is the free energy), the main control parameter is 
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P and the conjugate parameter of (3 with respect to —I3F is 

'-tm^-Eiff,j,Ly (2) 

The hnear series of equihbrium configurations in both ensembles, P{E) and —E{f3) at fixed J and L, 
are identical but the stability limits are not. Stability limits are obtained by applying the Poincare 
method to linear series (see Ledoux 1958). We shall consider the linear series of conjugate parameters 
(Katz 1978, 1979; Thompson 1979) since this is the most appropriate way to observe changes of 
stability. 

While state functions out of equilibrium, il, are never used, some assumptions must be made 
about J7 which are stated in section 2 where the extremization of the entropy is related to the Einstein 
equations when self-gravity is taken into account. Stability analysis in terms of Poincare's method is 
then briefly reviewed. We shall see that we need to know I3{E) only and not itself, a useful feature 
when analyzing the stability of self-bound radiation in a box. Furthermore, having determined fi{E) 
from the solution of the Einstein equations, we do not need to attribute a priori an entropy to the 
black hole. Instead, we recover it from the equilibrium configurations. 

The theory of fluctuations is then considered. Fluctuation theory as given by Landau and Lifshitz 
(1980) and Callen (1985) is not entirely applicable to self-gravitating systems for the reason that when 
gravitational interactions are important one cannot anymore split the system into a little subsystem 
and the rest which behaves like a heat reservoir for the little one. Nevertheless, we shall see that the 
mean quadratic fluctuations of the inverse temperature, in a small subsystem, are simply related to 
its own thermodynamic parameters as well as the parameters of the whole system. 

Section 4 gives details about (3{E) in non-rotating and rotating cavities. Fluctuations near turning 
points are calculated in section 5 and a summary of stability conditions and phase transitions in slowly 
evolving systems through a succession of quasi-equilibrium states is described in section 6. 

2. The stability of equilibrium configurations in mean field theory 

In a stationary axially symmetric distribution of matter in local thermodynamic equilibrium, with 
fixed total mass-energy and angular momentum, and with Einstein constraint equations given, the 
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total entropy S is stationary, 6S = 0, if Einstein "dynamical" equations are satisfied. Then Tolman's 
(1934) thermodynamic equilibrium conditions for local temperature hold and the angular velocity is 
uniform [Katz and Manor (1975); for non-rotating, spherical configurations see Cocke (1965)]. Here 
we shall use the partition function (Hawking 1978, Horwitz and Weil 1982, York 1988 and Brown 
et al. 1990) to sketch a formal deduction of both equilibrium and stability conditions for spherical 
distributions with zero angular momentum. This procedure will naturally exhibit the relations between 
the Einstein equations, quantum field theory in curved space and the extremization of the entropy. 

Consider, for definiteness, the state function Vt which encodes the total number of states of a scalar 
field in a curved spacetime with total mass-energy E as measured from infinity. The field is confined 
to a spherical box of "radius" L {i. e. the surface of the cavity is ^-nL?). Then il{E,L) is 

L) = e^(^-^) = Tr6{E - H) = j ^e^^* j |^P<^ exp [iSEin + iS^] ■ (3) 

The denominator ^^Diff" indicates that one should not integrate over geometries related by a diffeo- 
morphism; SEin = J d'^XyJ—gR, and 5*0 = / d^XyJ—gg^'^ 8^4)8^(1). Since there is no unique definition of 
a local matter Hamiltonian in general relativity, we have used the path integral formalism. We refer to 
Regge and Teitelboim (1974) for a definition of energy in asymptotically flat spacetimes. Integrating 
in (^) over all periodic matter configurations (of period t), one obtains 

L)= J ^J'^exp [iEt + iSE^n] Z[{9,,}- it] (4) 

where Z is the partition function at fixed Lorentzian time t in the background geometry {g^u}- 

Calculating (^) at the stationary configuration of lEt + iSsin + Z[{g^,y};it] gives the dominant 
contribution to 0,. At the stationary point, g^j^i, satisfies the time independent Einstein equations with 
an energy tensor of a thermal bath if at the saddle point t = —i(5. In the absence of horizons, the matter 
field configurations are defined for < r < L and one recovers, up to quantum matter corrections, the 
equations of a self-gravitating perfect fluid (Klein 1947, Gibbons and Hawking 1978). In the presence 
of horizons (see Car lip and Teitelboim 1993 for the treatment of the boundary term at the horizon), 
the regularity of the stress energy-momentum tensor, needed to satisfy Einstein equations, fixes the 
temperature of the matter (Hawking 1978). Then by integrating the time-time Einstein constraint 
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equation one obtains 



/ 



■r 



m{r, Mbh) = + 



dr'r'^{Tl{r'))HH,M.n 



(5) 



where {T^) HH,Mth the expectation value of the energy density, outside a black hole of mass M^h, in 
the Hartle-Hawking vacuum (Howard 1984). Since E = m{L,Mf,h), one may invert this relation and 
express M^h in terms of E and L. Therefore, one obtains the sought-for /? = (3{E) law (with L held 
fixed). Notice that Einstein equations with back-reaction taken into account provide an alternative 
way to determine the entropy of the black hole. Indeed, by integrating l3{E)dE one obtains the total 
entropy S{E,L), which reduces to the black hole entropy when the radiation energy is negligible. 
One may also subtract from S{E, L) the entropy of the radiation, but there is an arbitrariness in this 
subtraction — since there is no unique definition of the entropy density nearby the hole where the 
Tolman relation breaks down — which indicates that the concept of the entropy of an isolated black 
hole is probably meaningless. 

In order to illustrate the Poincare method, we perform the t-integration in (^) (Horwitz and Katz 
1978) which gives an 0, of the form 



We do not need to be specific about the function w{E, L, {g^^}) precisely because Poincare analysis 
deals with the succession of extrema of w{E, L, {Qfiu}) and not the function itself. A specific example 
is given in Sorkin et al. (1980) in the case of pure radiation without horizon, wherein w depends on 
Qrr only. 

The stationary value of w is the classical equilibrium entropy S. The local stability of equilibrium 
configurations is controlled by the quadratic fluctuations of w. Imagine we make a Fourier decompo- 
sition of the g^yS. Since the domain of existence is finite (0 < r < L), the g^^s are replaced by a 
denumerable set of variables, say, (i = 1, 2, ....). To order two, w may thus be written 



where an index e means 'at equilibrium' [thus S = We] and X''{E,L) are equilibrium values of x*. 
Let Aj's denote the (ordered: Ai < A2 < ...) elements of the diagonalized matrix —{d'^w/dx^dx^). 




(6) 




(7) 
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The Aj's are known as the Poincare coefficients of stabiHty (Lyttleton 1953). The equihbrium is 
thermodynamically (locally) stable, or stable for arbitrary small fluctuations, if and only if all the A^'s 
are positive i. e. if Ai > 0. In unstable situations, the number of Aj < characterize the degree of 
instability. 

The matrix —{d'^w/dx^dx^)e is a second order differential operator of dimension one over a finite 
domain. Wc may assume that the spectrum of eigenvalues is non-degenerate: Ai < A2 < •.• . Indeed, 
the slightest asymetric perturbation in a system would lift the degeneracy (Thompson and Hut 1973). 
Having then a non-degenerate spectrum of Poincare coefficients, the following properties hold (Katz 
1978, 1979): 

(a) Consider the linear series l3{E). Changes of stability along the linear series can only occur at ver- 
tical tangents, d/3/dE = ±00, like point A or C in figure 1. Such points where E is locally minimum 
or maximum are called turning points. Between two vertical tangents all equilibrium configurations 
have the same degree of instability. 

(b) In the neighborhood of a vertical tangent, when the linear series turns clockwise (its tangent goes 
from negative to positive values through infinity), one Poincare coefficient changes sign from nega- 
tive to positive value. That is, the system becomes stable or less unstable. If the linear series turns 
counter-clockwise, the changes of sign are reversed and the system becomes unstable or more unstable. 
It is thus enough to know the degree of stability of one configuration to know the degree of stability 
of all configurations. 

(c) Of particular interest are linear series with multiple turning points spiraling inwards against the 
clock. If we follow the spiral towards the limit point, we meet a succession of vertical tangents and 
beyond each one, an additional Poincare coefficient becomes negative. Counter-clockwise spirals rep- 
resent thus a succession of equilibrium configurations that are more and more unstable. 

Upon considering different ensembles, the following properties hold (Parentani 1994): 

(d) The most stable ensemble is always the most isolated one. The degree of stability of equilibrium 
configurations in any ensemble related by a Legendre transformation (which expresses the contact 
with a reservoir) to the most stable is immediately known if one knows the degree of stability in the 
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most stable ensemble. 



3. Fluctuations in gravitating systems 

In classical thermodynamics, it is well known (Landau and Lifshitz 1980 (§112), Callen 1985 and 
Landsberg 1990) that the mean square fluctuations of the fundamental thermodynamic quantities 
pertaining to any small part of a system (or to the system as a whole) are related to the specific heat. 



For instance, in a canonical ensemble of total volume V, the mean quadratic fluctuations (AE)^ of 
the energy E induced by the contact with the reservoir are given by 
d{-E) Cv 



= ^ = f (B) 



where Cy is the heat capacity at constant volume V . The mean square fluctuations (A/?)^ of the 
inverse temperature which result from those energy fluctuations are given by 



2 



since P is a function of E only. 

In a microcanonical ensemble, these fluctuations vanish since the total energy E is kept fixed. 
Nevertheless, within a small subsystem of volume V' , the temperature fluctuates and the mean square 
fluctuations of (3' are given by 



dp' /? 



12 



where Cy/ is the specific heat of the subsystem and E' its energy. This equation is valid only if the 
specific heat of the rest of the system is much bigger that the one of the subsystem (for homogeneous 
systems it requires V' <C V). The equivalence of ensembles (the equality of the fluctuations ( [lO| ) 
whether one works in the microcanonical or the canonical ensemble) means therefore that the rest of 
the system can be correctly treated as a heat reservoir for the little subsystem. Finally, we recall that 
the 'true' fluctuations are given by these estimates only if the characteristic dynamical time for the 
fluctuations to evolve is much bigger than f3 itself. 

In gravitating systems, there are long range forces. Therefore, when gravitational interactions are 
important, one cannot treat the rest of the system as a passive reservoir. Indeed the existence of stable 



microcanonical ensembles with negative specific heat indicates the decisive role of the energy constraint 
between the little subsystem and the rest. Furthermore, we stress the fact that when a microcanonical 
ensemble approaches instability, its heat capacity given in equation is always negative (thus the 
canonical ensemble is already unstable) for stable states and positive in unstable states (for which the 
canonical ensemble is still unstable) see section 4 and, for instance, figure 1. 

The relation between mean square fluctuations and thermodynamic coefficients in gravitating 
systems is thus more complicated than equation (^) and we shall display it in two steps. We shall 
first see that, since the vicinity of a turning point is dominated by a single eigenmode, one can relate 
the fluctuations of the least stable mode to the heat capacity of the whole system Cy- Then, we 
shall relate the fluctuations of the inverse temperature within a small subsystem to Cy itself. To prove 
the first point, consider a microcanonical ensemble (with J = 0) in which we define a temperature 
function (3{x'^;E,V) away from equilibrium (see also Okamoto et al. 1994) 

At equilibrium one has P{X^{E, V);E, V) = (3{E, V) where the equilibrium values x* = X''-(E, V) are 
solutions of 

(-) 

Therefore, the slope of the linear series and the derivatives of /3 are related by 

dp{E,V) _ f dP{x';E,V) \ sr-(dl\ dX'{E,V) 

dE \ dE i\dx'J^ dE ^ ' 

The derivative dX"^ /dE can be obtained from equation ( [l^ ) by the following identity 



d_ 
'dE 



dw 
dx^ 



' d~^\ dX' 



where there is no summation on i. Thus, from equation (|T4|) we deduce dX^ /dE in terms of {dp/dx'^)^ 
which we may insert into ( [T3| ) to obtain 

dp (dp\ {dl/dx/)l 

dE = [dE +2:^^- (15) 
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Consider now a linear series of stable configurations near a turning point. At that point the first 
eigenvalue Ai changes sign. As a result the right hand side of equation ( [T5| ) is entirely dominated by 
the first term 



(16) 



dE Ai \dx^ J ^' 

The changes of w are also dominated by the fluctuations of and given by 

u; ~ 5- ^Ai(Axi)2 (17) 



which, with equation ([T^), becomes 
2 dp/dE 



(18) 



where 

This A/3 represents fluctuations in the inverse-temperature function induced, near the critical point, 
by the fluctuations of the least stable mode x^. We can now use the standard arguments of fluctuation 
theory (Landau and Lifshitz 1980) and say that the probability dW for a fluctuation of /3 in the range 
P + A/3 and (3 + A/3 + d(3 is proportional to exp{w — S) and therefore, 

1 (A^)2 



(20) 



1 d(3 
V2^ VWJdE) ''''^ 2 {df3/dE) 

Thus the mean square fluctuations of /? are given by 

They are bounded when the microcanonical ensemble is stable (i. e. when Ai > or as we already 
point out when Cy < 0. See equation (|l^)). It is convenient to express dW in terms of Cy and of the 
dimensionless ratio 

P{x';E ,L) 

Then 



= u. (22) 



dW = J^e.v 



du (23) 



11 

and, following ( pl| ) the mean square fluctuations of u are 



{^uf = -Cy\ (24) 

Our task now is to relate those rather formal fluctuations of (3 to the fluctuations of the temperature 
within a small part of the system as well as to understand the origin of the unusual sign in equation 
(^l]). As we have already said, in the presence of long range forces, one cannot exactly split a system 
into two parts since there is no more a local definition of energy. If nevertheless some small part is 
less coupled gravitationally to the rest of the system, one may use it as the "small" system in which 
one can compute the fluctuations of 13. When this is not the case, we shall see that one can still 
consider the fluctuations within the outermost layer of radiation even though the energy into that 
layer is not well defined. We designate by the subscript 2 the little subsystem which is confined for 
radii L — I < r < L (where I <^ L). The subscript 1 refers to the rest of the system which is thus 
either the black hole surrounded by radiation up to that last layer, or pure radiation. Instead of the 
energy, we shall use the Schwarzschild mass (since it is a local quantity) and we thus split it into Ei 
and E — El. We shall see that the partition mass Ei will play a role very similar to the least stable 
mode x^. This is because d^iEi / 0. 

Let w{Ei; E, V) be the entropy out of equilibrium by which we now mean that the only variable 
out of equilibrium is Ei. All other variables have been replaced by their equilibrium values. Then 

w{Er,E,V) = Si{Ei,Vi) + S2{E-Ei,V-Vr,Ei) (25) 

where 5*1 and 52 are the entropies of the two parts. We emphasize the double dependence of Ei in 5*2. 
When El varies, the mass left over in the system 2 is E — Ei but the gravitational potential in which 
2 evolves is changed as well. This later dependence is nevertheless parametric if / ^ L. By parametric 
we mean that upon taking derivatives with respect to Ei, this later dependence gives an additional 
term which is 0{l/L) smaller than the usual term and which may be safely neglected. Equilibrium 
between the two parts requires, as usual, the equality of the temperatures: 

OeMEi-^E, V) = dE.Si + dE,S2 = /3i-p2 = (26) 

where Pi and P2 ^^e the inverse temperatures of the two parts and where we have neglected the 
additional term (This does not mean that we neglect completely the extra dependence in Ei, for 
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instance, ^2 depends explicitly on Ei trough the Tolman dependence (see the Appendix equation 
(A. 3))). Then the lowest order variations of w near equilibrium due to an exchange of energy AEi is: 



This energy fluctuation induces, in turn, a fluctuation of in the small part given by 



A/32 = ^E2{dE,P2) 



Then the entropy fluctuation induced by the latter one is 



(27) 



(28) 



w 



S = \{d(3E2) [l + {dE,P,){dpE2)\ {AP2? = l{dpE2)^idEMiAP2) 



(29) 



since the slope of the l3{E) curve of the entire system is given by 



dE 



{dE2p2 



dEi 
'dE 



(dE.P,) [1 + (dE^dp^^ 



Eo 



(30) 



because 



dE = dEi [l + {dp^E2){dE,Pi 



(31) 



Thus the mean square fluctuations of (32 given by the inverse coefficient appearing in equation (|2^ 
are 



d(3 



{A(32? = -{dEA)^{dpE^) 



(32) 



This is the generalisation of equation ( |10D when one cannot treat the rest of the system as a passive 
reservoir. Indeed, one does recover ( p!o| ) when the subsystem 1 can be approximated by the whole 
system in which case the coefficient of 9^2/^2 is 1- One sees also that the point wherein the fluctuations 
of (32 will diverge is controlled entirely by the divergence of d(3 /dE when the denominator in equation 
( |30|) vanishes. That is, one probes locally, through the fluctuations in the small subsystem the stability 
of the whole system 1 + 2. Furthermore, near the critical point, the mean fluctuations of /?2 are equal 
to the fluctuations of (i (and therefore independent of /) since (9^ E2){dEi(3i) = —1 at the critical 



point as seen in equation (^). This later factor of —1 explains the unusual sign in equation (21). 

For the interested reader, we also point out the analogy of equation (^) which gives the fluctuations 
when the rest of the system cannot be treated as a reservoir and the expression which relates the 
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fluctuations in one ensemble to the fluctuations in another ensemble related to the first one by a 
Legendre transformation (Parentani 1994). In both cases when the ensembles are nonequivalent, the 
fluctuations are controlled by the correction factor: -^{dp Ei), the coefficient of dE2P2 ™ equation 
(|32[). And it is only when the ensembles are equivalent that this factor reduces to 1. Another common 
feature is the fact that when d(3/dE = it does not imply that the fiuctuations of vanishes because 
dEiPi vanishes as well, see equation (pOl). 

For a canonical ensemble, when one consider configurations which are stable micro canonically, the 
necessary and sufficient condition for stability is the positivity of Cy- For those configurations only 
one may safely apply the analysis of Landau and Lifshitz and find that the fiuctuations of the total 
energy are indeed given by equation (^). Finally let us mention that equivalent ensembles would have 
vertical and horizontal slopes in (3{E) at practically the same point since both ensembles become 
simultaneously unstable. Thus the curve P{E) would make an angle bigger than 90°. Nothing like 
that is happening in the following situations. 



4. The linear series of conjugate parameters 

In Planck units, E, J and (3 are respectively given by 

E = ^, J=^, P = P*Ep (33) 
Ep n 

in which Ep = y/Kc^jG ~ lO^^erg and asterisks denote quantities in cgs units. Instead of using Planck 
units for E, J and /?, we shall use dimentionless quantities expressed in terms of G, c and L* as is 
common in classical general relativity. One thus define 

GE* _ GJ* _ nc(3* 

Then £ varies between and 1/2 (the Schwarzschild limit) and since J* < E*L*/c (the centrifugal 
limit), i7 is always smaller than £. Thus, 

0<J<£<^. (35) 
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We shall see nevertheless the appearance of a fourth dimensionless quantity L = L*/lp which will 
reflect the fact that the Hawking temperature has a quantum origin. Finally, the relation between E, 
J, 13 and £, J, b is 

(i) Schwarzschild Black Holes. 

Consider first a Schwarzschild black hole of mass MjJ'^ in equilibrium with radiation enclosed in a 
spherical cavity with fixed radius L. At sufficiently low temperature (high (3 or b), when most of the 
mass-energy is in the black hole, the total energy of the black hole and the radiation is, to a good 
approximation (see the Appendix) given by the sum of mass-energies as in flat spacetime (Hawking 
1976); in Planck units 

E = Ebh + Erad = ^ + ^nVp-^ (37) 
ovr 15 

where n is the sum over the helicity states of the massless fields. Translated into our classical units, 
( |37[ ) becomes 

47r^ / 1 \^ 

£ = £bh + £rad = b + cr = — — n ~ 6.91 • lO'^n. (38) 



45 \Stt J 

From now on we shall take n = 1. The function £{b) becomes eventful when £i,h and £rad are of the 
same order of magnitude. Indeed the lower bound for £ is reached when 

d£ , , 

at which point (C in figure la) one has 

£hh = ^£rad, £ = £c = \hc, bc = (4aL-2)i = 0.123 L"!. (40) 

When £ > £c-, the system's energy is dominated either by the black hole or by the radiation. 
Indeed for b > be, when £ > 1.6£c', less than 1% of the energy is in the radiation. For b < be, 
when £ > 25£c, less than 1% is in the black hole mass. At higher energies, along this later series, 
when b ~ L^^^^, one has to take into account the radiation self-gravitating effects and solve Einstein 
equations [Klein (1947) in the case of pure radiation, see also the Appendix]. One finds that there is a 
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maximal temperature where bo — O.llL"^/^ and Ed 2± 0.123 (see figure lb). For still higher energies, 
b increases again till one finds a maximal energy £a — 0.246 where 6^ ~ 0.135L~^/^ beyond which the 
are no more equilibrium configurations. 

Both curves, b{£) and b{£rad) are drawn in figure la and lb. In figure lb the two lines are 
almost on top of each other for £ ^ £r, (due to the small mass of the hole: £hhl^rad ~ 6 ~ L"^/^) 
and they spiral inwards counterclockwise with an almost common limit point Z at £z = 3/14 and 
bz ^ 1.236d ^ 0.132L-V2. 

Such counterclockwise spirals appear in Newtonian theory as well as in general relativity for equa- 
tions of state of the form P = Kp where P is the pressure, p the density and K a number (Chan- 
drasekhar 1972); here K = 1/3. The inward spiral will not differ very much whether we keep L or the 
proper volume of the sphere fixed. Therefore, stability limits will not differ significantly either. 

The two linear scries b[£) and b{£rad) coexist between two vertical lines for £c ^ £ ^ ^A- Thus 
between these lines there will be always one metastable state. We emphasis the origin of £c < £a 
which is due to the different nature of the instability: with or without a black hole at the Hawking 
temperature. This leads to the scaling 

|^ = 0.65L-2/^ (41) 

Very important also is the fact that there is a maximal temperature Td ~ 0.37L~^/^. Thus for 
L > 10^, T is always at least three orders of magnitude below the Planck temperature. We may 
thus safely neglect quantum matter effects (York 1985), quantum gravity and spontaneous nucleation 
(Piran and Wald 1984). For smaller boxes, not only should one take into account the quantum 
matter effects induced by the presence of the hot black hole, but one has presumably to abandon the 
thermodynamic analysis all together since the fluctuations become important and are not governed 
anymore by their thermodynamics estimates. Indeed, for those configurations, the thermal and the 
dynamical characteristic times come to coincide. 

Finally, we mention that Balbinot and Barletta (1988) computed the change of the black hole 
entropy S^h due to some quantum matter back-reaction. This provides a model for black hole remnants 
since the tiny black hole may now be in equilibrium with radiation for arbitrary large (or small) 
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temperatures. One easily shows that for cavities with L > 10^ one finds the hnear series to be 
almost unmodified and thus stability limits unaffected (contrary to what was suggested) since the 
modifications are dimensionalized by the Planck mass and since T ^ 0.37 10^^. 

(ii) Rotating Black Holes. 

Consider now a Kerr black hole at mid- height on the z axis of a cylindrical cavity of radius L*. 
When thermodynamical equilibrium is achieved, the radiation has the same temperature as the hole 
and rotates with the same angular velocity: = 0,rad- The peripheral velocity of the cavity, in units 
of c, is 

u = ^h^^ 0<u<l. (42) 



Schumacher et al. (1992) have calculated what becomes of equation ( pq ) when the gravitational pull 
of the hole and the self-gravitational effects of the radiation are neglected. This amounts to add to the 
energy of a rotating black hole £ijh{b,v) the energy of the radiation calculated (in special relativity) 
in rotating coordinates. We shall push the analysis of that approximation beyond its presumed limits 
of validity in order to gain some insight of what may still happen at very high angular momentum 
{v ^ 1) when general relativity has to be taken into account. 

Introducing a non-dimensional parameter h = J/Mbhrn where rjj is the radius of the horizon 
(Okamoto and Kaburaki 1992), the mass-energy, the angular momentum and the angular velocity of 
the hole can be written 

£th = Hl-''^), Jih = Y^{l-li^f-b\ 0<h<l (43) 

= (44) 

At equilibrium, r2f,/i = ^Icavj thus equations (42) and (^) give 



For a cylinder of height H* = 2L* , Schumacher et al. found 

. 1 - uys 

(1-1^2)2 



Sib, h) = 6(1 - h^) + |ctL-26-4 (46) 
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and 

J(5, h) = ^^{l- hYb' + aL-H-^j^^. (47) 



Equations ( ^7[ ) with (^) immediately indicate that most of the angular momentum is in the radiation 
except when b is close to 1/2, i. e. when the black hole fill up the entire cavity. For a given value of J 
held fixed, equation ( |47| ) gives h{b) which, substituted in equation (46) gives b{£). The b{£) function 
can only be written in parametric form and we have to solve it numerically. We have drawn b{£) for 
J' = 1/40 in figure 2a and for 1/8 in figure 2b. 

A common feature of the two linear series of figures 1 and 2 is that black holes and radiation coexist 
to the right hand side of a vertical line CCj-ad since there is always a minimum energy £c{-J) required 
to find a black hole in equilibrium with the radiation. There is no maximal energy, the equivalent 
of £a^ on the right hand side because we have not take into account the self-gravity of the rotating 
radiation. There are also two important modifications in the behavior of b{£) at fixed J when one 
compares linear series at different J . Both are shown in figure 3. First, the turning point C moves 
in the plane {£,b); £c increases with J. Second, the energy of the black hole at C, £bhciJ) — the 
energy gained by the radiation between C and Crad — divided by the total energy £c{J) decreases 
with J. 

Another novelty appears for J ^ 1/40. We remind the reader that in the absence of radiation, 
it was noted that for J^h/M^f^ ^ 0.68 the slope d/3/dMhh becomes negative, going through (Davies 
1977). When the radiation is taken into account, there is no change of sign in d(3/dE since the radia- 
tion dominates the equilibrium configuration before one approaches the critical ratio 0.68 of Jf,h/M^^. 
For J ^ 1/40 [] one does recover this phenomenon in a small interval of energy (see figure 2b). There 
are now two changes of sign upon decreasing the energy, before the black hole starts to evaporate. At 
lower energies the radiation dominates the equilibrium again and d(3/dE returns positive. 



5. Fluctuations near the turning points 

^ This limit is obtained by substituting v from (^5|) into (p7|). For a given J, one finds 3 values of < /i < 1 
corresponding to 3 points on h{£). Two of the points are on the left and right of the local maximum, point X in figure 
2b, when J is high enough, say 1/8. As J decreases the 2 points come closer to each other and point X goes down the 
line b{£bh)- At about J ~ 1/40 the points coincide and X is at the bottom of h{£bh)- For J ^ 1/40, there is no local 
maximum anymore. 
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We shall now compute the mean square fluctuations of the rescaled variable u, equation (|3^), near the 



turning points C and A. We shall also find the probabilities dW^ equation (23), that the fiuctuations 
be big enough for the system to jump from metastable to stable configurations. Finally we shall say 
a few words about characteristic times and for the probability rates. 

At the turning point Cy = 0. Thus, for an equilibrium configuration P in the vicinity of the 
turning point, to lowest order in [13 p — /3q) where (5q stands for either f3(j or one has 



where 



An = ^^^, |At.|«l. (49) 
Po 

Au parametrises the equilibrium configurations P and should not be confused with the fiuctuation 
out of equilibrium An defined in equation (p^). In terms of £ and 6, one has thus 



Cv = p-H%]^^-L-{h''^] An. (50) 



Near point C, equation (p8[) and (40) give 







- {Cv)c ^ l.gL^/^An. (51) 
Near point A, one finds 

- {Cv)a^23L^/\-Au) (52) 

(see the end of the Appendix for the numerical factor 23). The scaling law L^/^ itself is easily found: 
b{£rad) scales like L^^^"^ and the radiation entropy scales like cx L^^^. 

Near points A and C, the mean fluctuations of n, given in equation (p^, are extremely small (and 
therefore the temperature fluctuations as well) owing to the presence of the positive power of L in Cy- 
Nevertheless very close to the turning point, the fluctuations diverge and thus will be high enough to 
induce a transition towards stable conflgurations. We therefore estimate at which distance from the 
turning points, that is, for which Au, will the mean square fluctuations be high enough. To this end, 
we note that the minimal size that a fiuctuation should possess in order to provoke the transition is 
2Au. Indeed smaller fluctuations have not reached the minimal entropy configuration which lies onto 
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the unstable branch of the linear series and thus which is as well at a distance Au from the turning 
point. Hence the system will most likely return towards the initial metastable configuration. The 
minimal size of (Att)^ is thus reached when (Att)^ ~ (2An)^, that is when 

- Cy^ = A{Auf. (53) 

Using equation (|5^) we see that near point C this happens when {Au)c = 0.51L~^/^. At point A, 
following {Au)a ^ 0.22L~i/2. If L > 10^, one finds {Au)c < 2.0 • lO^^ and {Au)a < 2.2 • lO""^. 
Thus only very nearby the turning point will the mean fluctuations cause a phase transition. For Au's 
bigger than (Au)c or {Au)a the probabilities decrease drastically like exp(2(Ati)^Cy)). This shows 
that when Au is only a few times bigger than {Au)c or {Au)a the metastable black hole near C and 
the metastable radiation near A are perfectly stable. 

Having found dW, we now estimate probabilities of fluctuations per unit time. The rate at which 
a particular fluctuation of energy AE occurs is given, by virtue of the fluctuation dissipation theorem, 
by the inverse time it takes a cavity to return to equilibrium after the addition of the energy AE. This 
time depends on the peculiar dynamics of the system. For a black hole in equilibrium with radiation, 
the time, following Zurek (1980), is of the order of P'^AE. For pure radiation, the time is of the 
order of 0^ AE. Hence both times are proportional to powers of L. Thus probabilities per unit time 
behave essentially like dW itself, i. e. the negative exponential of powers of L dominates completely 
the probability rates. 



6. Stability and phase transitions in evolutionary sequences through quasi-equilibrium 
states: A summary 

(i) The Microcanonical Ensembles 

(a) Non-rotating Black Holes and Radiation. 

Consider an equilibrium configuration in the £-h plane, say at point F higher than point B in 
figure la. Let us remove energy by small amounts so that the system stays practically in equilibrium 
and evolution takes place along the linear series h = b{£). At energies £ < £b, the entropy of pure 
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radiation Srad is greater that the entropy of the composite system S. This is true only if one does not 
add a constant to the black hole Sjyh which will shift the point B. We emphasis that the addition of 
this constant, on the contrary, does not affect our computations of the mean fluctuations between B 
and C nor the fact that the black hole will evaporate at C. For £ < £b, the temperature is higher than 
the temperature at £b, the black hole is thus superheated (Gibbons and Perry 1978). We have just 
shown in section 5 that the probability for a fluctuation to lead to the total evaporation is completely 
negligible as long as Au ^ (Au)c — 0.51L~^/^. Thus the system will evolve almost down to £c 
staying in these superheated states. If one remove energy below £c, the black hole cannot survive in 
equilibrium anymore and will evaporate into radiation to a 6 = bcrad = 0.082L~^/^, that is, with a 
temperature given by B^^'^Tc, since 4/5 of the energy was in the black hole. Removing more energy 
will simply cool the radiation down. 

Now if one reverses the process and starts to add energy to the cavity, say from point G in figure 
la, the evolution takes place along the linear series b{£rad)- The radiation heats up and reaches the 
point Brad where £Brad = ^B- Beyond this point the radiation finds itself also in a superheated state. 
The chances to form a black hole at those energies are exponentially small (see also Piran and Wald 
1982). The radiation will thus continue to evolve (from figure la to figure lb) in that superheated state 
as long as \A.u\ ^ |Au|^ ~ 0.22L~^/^, that is, almost up to point A. It will never become supercooled 
because 6^ ^ 0.14 L~^/^ < hsrad — 0.20 for L > 10^. Near point A, the radiation will collapse 
and form, near point A^h, a black hole in equilibrium with the left over radiation. 

One has thus a closed circuit of equilibrium configurations which can be experienced counterclock- 
wisely only. 

(b) Black-Holes in Rotating Cavities — Self-Gravity Neglected. 

If we compare figures 1 and 2, we see that rotation does not modify very much the above picture. 
Even for fast rotation, significant differences occur only in the late stages of evolution before evapo- 
ration. As one approaches point C in figure 2, the black hole has already lost most of its energy (see 
figure 3) . The higher the angular momentum, the smaller this energy. In evolutionary sequences with 
decreasing energy of fast rotating cavities, evaporation of black holes will go almost unnoticed, the 
first order phase transition being very mild. Since we have not taken into account the self-gravity of 
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the radiation, one does not find the equivalent of the point A nor the closed circuit. 

(ii) Canonical Ensembles 

(a) Non-rotating Black Hole in a Cavity 

If one assume that one may control the temperature at infinity instead of letting the system be 
isolated, stability conditions can be read from the same b{£) diagram but rotated 90° clockwise (or 
equivalently one looks for horizontal tangents). Looking at figure 1, we observe the following situation. 
At low temperature, high b, with very little radiation, the lonely black hole cannot be stable since it has 
negative heat capacity. Then the —S{b) curve rotates counterclockwise and the sequence of inwards 
spiraling configurations will become more unstable each time one encounters a vertical tangent. Thus 
black holes in a cavity with radiation at fixed b are always unstable. We mention that York (1985) has 
considered cavities with fixed temperature at the boundary rather than fixed temperature as measured 
from infinity (see equation (A.3) in the Appendix). This leads to a different Legendre transformation 
which defines another free energy. Therefore, it is not surprising that this other ensemble has different 
stability limits when the energy approaches the Schwarzschild limit. 

One should, however, question the physical relevance of canonical configurations in general wherein 
self-gravitating effects are important. This is because one cannot ignore the self-gravity of the reservoir 
needed to fix the temperature. Indeed, in order to maintain properly the temperature fixed, the 
reservoir has to be large compared to the system itself. Thus the reservoir would be within its own 
Schwarzschild radius. It appears therefore that the canonical situations have hardly any physical 
relevance. 

Pure radiation behaves quite differently. Black body radiation in a cavity at low temperature (and 
hence small mass) is stable. If we slowly heat the cavity, the sequence of evolution is the —£rad{b) 
curve of figure 1. Equilibrium configurations are all stable up to a temperature Td or down to point 
D where Cy becomes negative. For b below bo there is no equilibrium configuration and at higher 
energies the configurations are all unstable, since —£rad{b) spirals inwards counterclockwisely. 

(b) Black Holes in Rotating Cavities — Self-Gravity Neglected 

Black holes in rotating cavities and in a heat bath are as dull as non-rotating black holes. They 
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are all unstable. The situation changes, however, somewhat when J ^ 1/40. Then, as can be seen 
on figure 2b rotated clockwise 90°, two vertical tangents appear at point X and Y enclosing a narrow 
range of energies. At high 6, i. e. b > bx, the canonical ensemble is certainly unstable, the specific heat 
being negative. Since we know that the micro canonical ensemble is stable for those configurations, the 
number of negative Poincare coefficients is only one (see point (d) of section 2). From the point X, the 
system will stay stable if we slowly change the temperature as to decreasing the energy from £x to 
£y- At point Y we reach again a vertical tangent. The counterclockwise turn of the line signals that 
instability is back. All equilibrium configurations at higher temperatures are thus unstable. Notice 
that the "stabilizing" effect of the angular momentum appears in a domain where general relativity 
begins to be important. 
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Appendix A: On self-gravitating radiation in a box 

The properties of global equilibrium configurations of a self-gravitating, finite-sized, non-rotating, 
spherical symmetric fluid with no black hole in the center has been studied for various reasons by 
Klein (1947), Chandrasekhar (1972), Sorkin et al. (1982) and Page (1992). The equation of state of a 
gas of photons (n = 1) is 



where p is the mass-energy density. Under local thermodynamic equilibrium conditions [see equation 
(|8|)], p is given by 




(Al) 




{A.2) 



in Planck units where T(r) is the local temperature function of r, T{r) = TgoQ{r) where is the 
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r-dependent lapse function which relates the local proper time to the asymptotic one. In particular, 
on the surface of the cavity, one has 



r(L)^l - "^^radiL) ^ TiL)^l-2£rad{L) =T=^. (A3) 

where T is the "temperature as measured at infinity" . The distribution p{r) is given by the Tolman- 

Oppenheimer-Volkoff equation which can be solved as follows (Sorkin et al. 1981). Set 

q = ATvr^p, l^ = ¥ jl dr'r'^p'. {AA) 

Then the TOV equation reads 

dq _ 2g(l-4//-(2/3)g) . 
(l-2/x)(g-/x) • ^^-^^ 

Initial conditions are ^ = q = 0, this leads to a family of solutions parametrized by the density at 

r = 0. On the boundary r = L 

q{L) = AirL^piL), p{L) = S^ad- (^-6) 

The second equation is used to express the density at r = in terms of £rad- Thus, the solution of 

equation (A. 5) defines, with (A. 3) and (A. 6) the function b{£rad) 

q{l - 2£radf ^ ' 

where g is a function of £raAi- 

When self gravity is negligible ~ 3/i ^ 1), h{£rad) reduces to the expression valid in flat space 
time E^aci = c^VT^. When £rad > 0.1, self gravity is important and £rad{b) is given by (A.8). Special 
points are given in the table (A.9); the limit point Z of figure 1 is where the density at r = is infinite. 



Special points* 


D 


A 


E 


Z 


£rad — 


0.123 


0.246 


0.235 


0.214 




0.107 


0.135 


0.140 


0.132 




0.372 


0.295 


0.284 


0.301 



(A.9) 



'See figure 1. For comparison, £c 0.154 L and bc^^^ ^ 0.082L 



When a black hole is present in the center of the cavity, (A.l) and (A.2) are no more valid. 
One has to take into account the quantum behavior of the mass energy density since the Tolman 
relations certainly break down near the horizon. The quantum version of the density is provided by 
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the expectation value of the time-time component of the stress energy tensor in the so-called Hartle- 
Hawking "vacuum" (Howard 1984). One finds first that the contribution to the gravitational mass for 
2Mbh <r < QMhh is of the order of M^"^^ ( in our units, it means b ) and secondly that for r > GM^/j 
one may approximate the mass energy density by (A. 2). Thus one may approximate (|5|) by 

m{r, Mth) = Mkh + Qj,, dr'r'^M^^{l - 2M^,h/r')-\ (.4.10) 

(again for cavities with L > 10^). We have also neglected in (A. 10) the self-gravity of the radiation. 

For b > be, this is certainly valid since Erad < E/5 and since Eq <^ L/2. For smaller black holes 
{Efjfi < (4:/5)E), the configurations are unstable and it is therefore useless to calculate the corrections. 
From (A. 10) and for b > be, it is easy to verify that the corrections to b{£) are small and have no 
effect on thermodynamic stability limits. The reader might consult Page (1992) to find an explicit 
evaluation of corrections to the entropy. 

We now give the derivation of {Cv)a appearing in equation ([5^). We use equation (A. 8) with 



(All) 



equation (A. 6) and obtain the following expression for d£/db (we now drop the index rad) 
3aL-U£ ^ q{l-2£f{q-£) 

2{^-g-l + 2£) ■ 

Using again (A. 6) with (^9]), we also obtain for dq/db 

2 . 



q\l-2£y{l-i£ 



A. 



SaLZ^dq _ 

(|,-l + 2f) 

These two equations give us a means to calculate the second derivative of d'^S/db'^ 



(A12) 



(SaL-^d'^q 
\ 46^ db' 



d 
'Uq 



q{l-2Sf{q-£) 
2{^q-l + 28) 



d6 



{A.13) 



At point A, the first derivative is zero 
{-Cv)a = 64iL2(i^£lfl4^Aa, 

With bji given in the above table, one has thus 
(-Cy)A^23L3/2(-Au). 



(^.14) 



(A15) 
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Figure captions 

Figure 1. b{6bh), b{6rad) ^'^d b{S) are drawn for L = 10^, in figure la for S ^ 10~^ and in figure lb 
for S ^ 10-^ Between S = 10'^ and 10-^ b{Srad) and b{S) come closer and closer to the point that in 
figure lb the two lines are indistinguishable. The line b{£) through points FCQDAZ is a counterclock 
inward spiral. The dotted line in the lower left hand corner of figure lb is the non-relativistic (no 
selfgravity) b{Srad) for comparison. 

Figure 2. b{£bh,J'), b{£rad-,J) and b{£,J) curves for constant angular momentum J and L = 
2.65 • 10^. The low value of L is to make clear figures. Figure 2a is for J = 1/40 and figure 2b is 
for J = 1/8. Both lines are drawn with the same limits of h and £, showing the displacement to 
right for increasing J of the linear series b{£rad-> J)-> b{£, J7) and of the point C. Once the equilibrium 
configurations leave the linear series b{£bh,J'), the black hole starts to loose mass with respect to the 
radiation, see figure 3. 

Figure 3. This figure displays two lines: (a) The minimum of energy £c{J) as a function of increasing 
angular momentum. The scale of £c displayed on the left is the same as in figure 2. The curve is 
parametrized in real values of £c- Notice that beyond J = 3/16 ~ 0.2, the energy £c is highly 
relativistic and the non-relativistic curve is likely to be different, (b) The ratio £bhc/£c = bc/£c 
with its scale displayed on the right. For J" = 0, £\,hcl^C = 4/5; this is well above the limits of the 
drawing. 
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